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Abstract 



In this paper, we characterize the odd positive integers n satisfying 
the congruence Y^=\ 3 _5_ — (mod n) . We show that the set of 
such positive integers has an asymptotic density which turns out to 
be slightly larger than 3/8. 



1 Introduction 

Given any property P satisfied by the primes, it is natural to consider the 
set Cp := {n composite : n satisfies P}. Elements of Cp can be thought of 
as pseudoprimes with respect to the property P. Such sets of pseudoprimes 
have been of interest to number theorists. 

Putting aside practical primality tests such as Fermat, Euler, Euler- Jacobi, 
Miller-Rabin, Solovay-Strassen, and others, let us have a look at some in- 
teresting, although not very efficient, primality tests as summarized in the 
table below. 





Test 


Pseudoprimes 


Infinitely many 


I 


(n — 1)! = —1 (mod n) 


None 


No 


2 


a n = a (mod n) for all a 


Carmichael numbers 


Yes 


3 


YT^3 m = -1 (mod n) 


Giuga numbers 


Unknown 


4 


<j>{n)\{n-l) 


Lehmer numbers 


No example known 


5 


E"=ii n_1 = - 1 ( mod n ) 




No example known 



In the above table, <f>(n) is the Euler function of n. 

The first test in the table, due to Wilson and published by Waring in |T9] . 
is an interesting and impractical characterization of a prime number. As a 
consequence, no pseudoprimes for this test exist. 

The pseudoprimes for the second test in the table are called Carmichael 
numbers. They were characterized by Korselt in [10]. In [I], it is proved that 
there are infinitely many of them. The counting function for the Carmichael 
numbers was studied by Erdos in [6] and by Harman in [9]. 
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The pseudoprimes for the third test are called Giuga numbers. The sequence 
of such numbers is sequence A007850 in OEIS. These numbers were intro- 
duced and characterized in [3]. For example, a Giuga number is a squarefree 
composite integer n such that p divides n/p — 1 for all prime factors p of 
n. All known Giuga numbers are even. If an odd Giuga number exists, it 
must be the product of at least 14 primes. The Giuga numbers also satisfy 
the congruence nB^ n ^ = — 1 (mod n), where for a positive integer m the 
notation B m stands for the mth Bernoulli number. 

The fourth test in the table is due to Lehmer (see [IT]) and it dates back to 
1932. Although it has recently drawn much attention, it is still not known 
whether any pseudoprimes at all exist for this test or not. In a series of 
papers (see [13], [15], and [IS]), Pomerance has obtained upper bounds for 
the counting function of the Lehmer numbers, which are the pseudoprimes for 
this test. In his third paper [IB], he succeded in showing that the counting 
function of the Lehmer numbers n < x is 0(a; 1/ ' 2 (logx) 3//4 ). Refinements 
of the underlying method of [16J led to subsequent improvements in the 
exponent of the logarithm in the above bound by Shan [T7j, Banks and 
Luca [2], Banks, Giiloglu and Nevans [3], and Luca and Pomerance [12], 
respectively. The best exponent to date is due to Luca and Pomerance [T2] 
and it is —1/2 + e for any e > 0. 

The last test in the table is based on a conjecture formulated in 1959 by 
Giuga [S], which states that the set of pseudoprimes for this test is empty. 
In 0|, it is shown that every counterexample to Giuga's conjecture is both a 
Carmichael number and a Giuga number. Luca, Pomerance and Shparlinski 
[T3] have showed that the counting function for these numbers n < x is 
0(x 1//2 /(logx) 2 ) improving slightly on a previous result by Tipu [18J. 

In this paper, inspired by Giuga's conjecture, we study the odd positive 
integers n satisfying the congruence 



It is easy to see that if n is an odd prime, then n satisfies the above congru- 
ence. We characterize such positive integers n and show that they have an 
asymptotic density which turns out to be slightly larger than 3/8. 




(mod n). 



(1) 
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For simplicity we put 

n-l 

gw = £j L( "- 1)/2J , 

although we study this function only for odd values of n. 



2 On the congruence G(n) = (mod n) for odd 
n 

We put 

:= {n odd : G(n) = (mod n)}. 
It is easy to observe that every odd prime lies in *JJ. In fact, by Euler's 
criterion, if p is an odd prime, then j( p_1 )/ 2 = ( - ] (mod p), where ( - J 
denotes the Legendre symbol of j with respect to p. Thus, 

p-i 



so that p e 

We start by showing that numbers which are congruent to 3 (mod 4) are in 

qj. 

Proposition 1. Ifn = 3 (mod 4), £/ien n G *}5. 

Proof. Writing n = 4m + 3, we have that (n — l)/2 = 2m + 1 is odd. Now, 

n-l 

2G(n) = ^(j 2m+1 + (n-j) 2m+1 ) 

n-l 

= n {f m + 3 2m - X (n-j) + --- + (n- j) 2m ) , 

so n | 2G(n). Since n is odd, we get that G(n) = (mod n), which is what 
we wanted. □ 
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The next lemma is immediate. 

Lemma 2. Let p be an odd prime and let k > 1 be an integer. Then 



p k — 1 , k\\ i f P k ~ 1 iA J p — 1 if k is even, 

(p — l)/2 «/ zs ooW. 



gcd — - — ,<p(p ) = gcd — - — ,p - 1 



With this lemma in mind we can prove the following result. 

Proposition 3. Let p be an odd prime and let k > 1 be any integer. Then, 
p k G ^3 if and only if k is odd. 

Proof. Let a £ Z be an integer whose class modulo p fc is a generator of the 
unit group of Z/p fc Z. We put (3 := a^ p ~ 1 ^ 2 . Suppose first that k is odd. We 
then claim that /3 — 1 is not zero modulo p. In fact, if a^ p ~ 1 ^ 2 = 1 (mod p), 
then since also a p ~ l = 1 (mod p), we get, by Lemma [21 that a^ 1 ^ 2 = 1 
(mod p), which is impossible. 

Now, since — 1 is coprime to p, it is invertible modulo p k . Moreover, since 
also k < (p k — l)/2, we have that 

n-1 

Gf(n) = Y,J ipk ' 1)/2 = Yl 3 ipk ~ 1)/2 ( mod ^) 

3=1 gcd(i,p)=l 
l<i<n-l 



^ (aW-*)/*y (mod p *) = ^ ? (mod p fc ) 
(modp fc ). 



j=l i=l 



0-1 

Assume now that is even. Observe that 

(p fc - l)/2 = (p - 1)((1 + p + • ■ • + p fc ^)/2) := (p - l)m, 

and m is an integer which is coprime to p. Thus, (3 = -1 )/ 2 = (a( p_1 )) m 
has order p fc_1 modulo p fe , and so does a v ~ x . Moreover, again since k < 
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(p — l)/2, we may eliminate the multiples of p from the sum defining G(n) 
modulo n and get 

G(n) = J2j(p k -m= J {pk - 1)/2 (mod/) 

3=1 gcd(j,p)=l 
l<j'<n-l 

<p(p k ) ■ p fe_1 0-i) 

= J2 (a^'^y = (« {p_1) ) im (mod/) 

i=i ' i=i 

p k-i 

= (p-l)^(^- 1 ) 4 (mod/). (2) 

i=l 

Since has order p k ~ 1 modulo p h , it follows that a p ~ x = 1 + pu for some 
integer u which is coprime to p. Then 

£V')' = «(^)- (3) 

Since a pfe_1 = 1 +/« (mod / +1 ), it follows that (a pk ' 1 — l)/(a — 1) = p k ~ l 
(mod p k ), so that 

= ap fc_1 (mod/)=p fe-1 (mod/). (4) 

Calculations (j3J) and (jlj) together with congruences (jSJ) give that G(n) = 
(p — l)p (mod p ). Thus, p is not in ^ when is even. □ 

Note that Proposition [3] does not extend to powers of positive integers having 
at least two distinct prime factors. For example, n = 2021 = 43 x 47 has the 
property that both n and n 2 belong 




3 A characterization of *p and applications 

Here, we take a look into the arithmetic structure of the elements lying in 
^3. We start with an easy but useful lemma. 
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Lemma 4. Let n = Y[ p r p\\nP rv be an odd integer, and let A be any positive 
integer. If gcd(A,p — 1) < p — 1 for all p \ n, then 

j A = (mod n). 

gcd(j»=l 
l<i<n-l 

Proof. It suffices to prove that the above congruence holds for all prime 
powers p rp ||n. So, let p r be such a prime power and let a be an integer which 
is a generator of the unit group of Z/p r Z. Put (3 := a A . An argument similar 
to the one used in the proof of Proposition [3] (the case when k is odd) shows 
that the condition gcd(A,p — 1) < p — 1 entails that — 1 is not a multiple 
of p. Thus, f3 — 1 is invert ible modulo p. We now have 

0(p r ) 



£ / = (^r) £ ^ (mod p-) ^ #n/jO £ a* (mod p 

/„• „\_i VW // g C( jc„- , i 

1; 

<Kp r ) .,6iV)+i _ o 



gcd(j,n)=l gcd(j,p)=l i=l 

l<j<n-l l<i<P 



= (f)(n/p r ) f3 l (mod p r ) = (f)(n/p r ) (mod p r ) 

= (mod p r ), 

which is what we wanted to prove. □ 

Theorem 5. A positive integer n is inty if and only ifgcd((n — l)/2,p— 1) < 
p — 1 /or a// p | n. 

Proof. Assume that n is odd and gcd((n — l)/2,p — 1) < p — 1. By Lemma 
H 

£ j (n ~ 1)/2 = (modn). 

l<7<n-l 

Now, let ei be any divisor of n. Observe that 

J 2 = d 2 ^ Z 2 . (5) 

(j,n)=d (i,n/d)=l 
l<i<n-l l<i<n/d-l 

The last sum in the right-hand side of fl5]) above is, by Lemma HI a multiple 
of n/d, so that the sum in the left-hand side of §5§ above is a multiple of n. 
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Summing up these congruences over all possible divisors d of n and noting 
that 

G{n) = Yl E ^ 1)/2 > 

d\n gcd(j,n)=d 
l<j<n-l 

we get that G{n) = (mod n), so n G 

Conversely, say n G ^ is some odd number and assume that there exists a 
prime factor p of n such that p — 1 | (n — l)/2. Write (n — l)/2 = (p — l)m. 
Observe that m is coprime to p. Assume that p r \\n. Then, modulo p r , we 
have 

n-1 

G(n) = £ j^/ 2 = (n/jf) E 3 {n ~ 1)l2 (mod /) = (n/pO £ j^. 

i=l gcd(j,p)=l gcd(j,p)=l 

l<j<P r -l l<J<P r -l 

The argument used in Proposition [3] (the case when k is even) , shows that 
the second sum is not zero modulo p r , and since n/p r is also coprime to p, 
we get that p r does not divide G(n), a contradiction. 

This completes the proof of the theorem. □ 
Here are a few immediate corollaries of Theorem [5j 

Corollary 6. Let n be any integer. Assume that one of the following condi- 
tions hold: 



i) gcd ((n — l)/2, ip(n)) is odd; 

ii) gcd ((n — l)/2, A(n)) is odd, where A(n) £/ie Carmichael function. 
Then n G 

Corollary 7. 7/n fe G ^3 /or some > 1, i/ien nG^. 

Proof. Observe that gcd ((n — l)/2,p — 1) divides gcd ((n fc — l)/2,p — l) for 
every and every prime number p. Now the corollary follows from Theorem 
El □ 
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We add another sufficient condition which is somewhat reminiscent of the 
characterization of the Giuga numbers. 

Proposition 8. Let n = Yl p r p ^ n p rp be an odd integer. Ifp — 1 does not divide 
n/p rp — 1 for every prime factor p of n, then n G 

Proof. By Theorem [5j if n then there exists a prime factor p of n such 
that p — 1 divides (n — l)/2. In particular, p — 1 | n — 1. Since p — 1 also 
divides p rp — 1, it follows that p — 1 divides n — p rp = p rp (n/p rp — 1). Since 
p — 1 is obviously coprime to p rp , we get that p — 1 divides n/p Vp — 1, which 
is a contradiction. □ 

It is also easy to determine whether numbers of the form 2 m + 1 are in 
Indeed, assume that 2 m + 1 G" for some positive integer m. Then, by 
Theorem El there is some prime p \ 2 m + 1 such that p — 1 \ ( (2 m + 1 ) — 1 ) /2 = 
2 m ~ 1 . Thus, p = 2 a + 1 for some a < m — 1, and so p is a Fermat prime. In 
particular, a = 2 a for some a > 0. Since p = 2 2 " + 1 is a proper divisor of 
2 m + 1, it follows that 2 a | m and m/2 a is odd. This is possible only when 
2 a is the exact power of 2 in m and m is not a power of 2. So, we have the 
following result. 

Proposition 9. Let n = 2 m + 1 and m = 2 a m\ with a > and odd mi > 1. 

Then n 6 unless 2 2 ° + 1 is a Fermat prime. 



4 Asymptotic density of 

Let I be the set of odd positive integers. In order to compute the asymptotic 
density of ^)3, or to even prove that it exists, it suffices to understand the 
elements in its complement It turns out that this is easy. For an odd 

prime p let 

JF p := {p 2 (mod 2p(p — 1))}. 

Observe that T p C I. 
Theorem 10. We have 

IMP = |J T o- (6) 

p>3 
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Proof. By Theorem [5J we have that n if and only if p — 1 divides 

(n — l)/2 for some prime factor p of n. This condition is equivalent to n = 1 
(mod 2(p — 1)). Write n = pm for some positive integer m. Since p is 
invertible modulo 2(p — 1), it follows that m is uniquely determined modulo 
2(p — 1). It suffices to notice that the class of m modulo 2(p — 1) is in fact 
p since then pm = p 2 = 1 (mod 2(p — 1)) with the last congruence following 
because p 2 — 1 = (p — l)(p + 1) is a multiple of 2(p — 1). This completes the 
proof. □ 

Observe that J-p is an arithmetic progression of difference l/(2p(p—l)). Since 
the series 



is convergent, it follows immediately that I\^P; hence, also has a density. 
This also suggests a way to compute the density of *}3 with arbitrary precision. 
Namely, say e > is given. Let 3 = p\ < p 2 < ■ ■ ■ be the increasing sequence 
of all the odd primes. Let k := k(e) be minimal such that 



It then follows that numbers n ^ ^ which are divisible by a prime pj with 
j > k belong to {Jj >k J r Pj , which is a set of density < e. Thus, with an error 
of at most e, the density of the set I\*}3 is the same as the density of 





j<k 



which is, by the Principle of Inclusion and Exclusion, 






11,12, ■■;is 



(7) 



s>l l<ii<i2< - <«s<fe— 1 



with the coefficient e^,^,...,^ being zero if f] s t=1 J r p i = 0, and being (— 1) 
otherwise. Taking e := 0.00082, we get that k = 29, 



274510632303283394907222287246970994037 
2284268907516688397400621108446881752020 



0.120174, 



j<29 
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and consequently belongs to [0.379005, 0.379826]. So, we can say that 

p(*P) = 0.379 . . . 

Here and in what follows, for a subset A of the set of positive integers we 
used p(A) for its density when it exists. 

These computations were carried out with Mathematica, for which it was 
necessary to have a good criterion to determine when the intersection of T v 
for various odd primes p is empty. We devote a few words on this issue. Let 
us observe first that the condition n 6 T v , which is equivalent to the fact that 
p | n and p — 1 divides (n — l)/2, can be formulated as the pair congruences 

n = 1 (mod2(p— 1)); 

n = (modp). (8) 

Assume now that T 3 is some finite set of primes. Let us look at f] p€V T v . Put 
m := ripG-pP- The ^ rs ^ se ^ °f congruences (|8]) for all p EV is equivalent to 

n = 1 (mod 2A(m)), (9) 

where A(m) = lcm [p — 1 : p G "P] is the Carmichael A-function of m. The 
second set of congruences for p E V is equivalent to 

n = (mod m). (10) 

Since 1 is not congruent to modulo any prime q, it follows that a necessary 
condition for (jUJ) and (JTUj) to hold simultaneously is that m and 2A(m) are 
coprime. This is also sufficient by the Chinese Remainder Lemma in order 
for the pair of congruences (JSJ and ( llOp to have a solution n. Since m is also 
squarefree, the condition that m > 1 is odd and m and 2A(m) are coprime 
is equivalent to m > 2 and m and 0(m) are coprime. Put 

M := {m > 2 : gcd(m, 0(m)) = 1}. (11) 

Thus, we proved the following result. 

Proposition 11. LetV be a finite set of primes and put m := ripe^p^- Then 
Op^pFp is nonempty if and only if m 6 Ai, where this set is defined at ( ITT]) 
above. If this is the case, then the set f] pe p T v is an arithmetic progression 
of difference 1 / (2mA (m) ) . 
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The condition that m <E M can also be formulated by saying that m is odd, 
squarefree and p \ q — 1 for all primes p and q dividing m. We recall that the 
set Ai has been studied intensively in the literature. For example, putting 
M(x) — M fl [1, x], Erdos [5] proved that 

#M(x) = e~ 7 (l + o(l))- — as x^oo. 

log log log X 

In particular, it follows that if V is a finite set of primes, then f] peV T v ^§ 
if and only if T v f] J- q ^ for any two elements p and qoiV. 

Finally, let us observe that with this formalism and the Principle of Inclusion 
and Exclusion, as in ([7]) for example, we can write that 

PW) = 1^ 2m\{m) ' 

m<=MU{l} y ' 

Here, u(m) is the number of distinct prime factors of m. The fact that the 
above series converges absolutely follows easily from the inequality X(m) > 
(logm) clogloglogm which holds with some positive constant c for all sufficiently 
large m (see [7]), as well the fact that the series 

E 



m(logm) 2 

■m>2 v to ' 



converges. We give no further details. 
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